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Abstract 
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1 Introduction 



The status of the theoretical description of the Drell-Yan-hke processes p] was widely overviewed in the 
resent papers [2j [3] where the necessity of further in-depth study of them both from experimental and 
theoretical sides was emphasized. 

Single W and Z boson production in Drell-Yan-like (DY) processes have clean signature and large 
cross sections. Their precision studies will be used for determination of the parton density functions, 
luminosity, M^, sin^ 9^^^^ , ^w,z- Required precision tag is below 1%, see e.g. [HIS]. 

Theoretical calculations of the DY processes for high energy hadronic colliders were performed at the 
level of one-loop QED and electroweak (EW) radiative corrections (RC) by several groups, see papers 
[71 [HI [HI [ini [H] and references therein. QCD corrections are known up to the next-to-next-to-leading 
order [T^ . 

The first attempt to combine EW and QCD corrections was done within the working group "The 
neutral current Drell-Yan process in the high invariant mass region" of the the 2007 Les Houches Work- 
shop [13]. Existing codes have to meet the requested experimental precision. QED, Weak and QCD 
one-loop level corrections should be calculated taking into account by oneself, their interplay and all 
necessary higher order effects. 

For users attention we offer two types of SANC outputs 0: 

• stand-alone packages for calculation of the EW and QCD NLO RC at the parton level together with 
the environment in which it is run, i.e. the Standard SANC FORM (and/or FORTRAN) Modules 
(SSFM) (see section [SHI; 

• MC event generators, for production of event distributions at the hadronic level, based on the FOAM 
algorithm [TO] . 

In this note we present the first positive experience of the using SSFM's for DY-like processes for 
neutral current (NC), 

qq ^ l+e- (1) 

with £ = ij,,e, (the NC DY generator was presented first in talk [XT') and charged current (CC) processes 
(in — p]) 

qq' ^ ei^e , (2) 

as stand-alone codes, see also [19] and [20] . 

The packages for calculations at the partonic level as well as the FOAM-based generators are accessible 
from SANC project homepages [21]. 

Moreover, we create tools for checking the implementation of these modules — the integrators of the 
processes, based on the Vegas algorithm l22]-[23], but presently only for an internal use. We used these 
integrators in the business of the tuned comparison within the several international workshops for NC 
and CC DY processes, see Les Houches Workshop Proceedings [23], [13] and Tevatron for LHC Report 

m- 

In this note we demonstrate the first example of applying EW SSFM into MC generator WINHAC 26J. 
We checked the work of these modules by means of comparison of the results obtained by SANC integrator 
DY_CC_VEGAS with those produced by the development version WINHAC vl.30 (see Section[3]). 

The SANC project was presented at several ATLAS MC and SM working group meetings, e.g. [?7l [251 
[221 [3ni [35 • From the project page one can download SANC client software, that allow users to interactively 
test the workability of the all offered processes, i.e. user may create SSFM's by oneself (the documentation 
for SANC client also can be found at the project webpage). However, user may download SSFM from the 
SANC homepages directly. 

^We do not describe SANC system in this note, referring the reader to pubhshed papers |14| . |15| . 
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2 Basics of working with SANC 



2.1 Working with SANC client 



To work with SANC, one must install the SANC client. The SANC chant can be downloaded from the SANC 
project homepage [H]. On the homepage select Download, see FigH] then download the client, unzip 
it and follow the instructions in the README filell 



SUPPORT OF ANALYTIC AND NUMERIC CALCULATIONS FOR EXPERIMENTS AT COLLIDERS 

New Clientv.l.lO ay/aitable! 



SANC Project 



I ^iews 



Team Members | Users SjppDrt~ 



SANC Project Download. 



:^ SANC Client -- Speciat interface for networking at SArfC Server. 

±' SANC Modules -- Stand aione SANC products for precision cakulations. 

*' SANC Generators -- Stand alone SANC packages for events generatfon. 



SANC Client. 

Combined 5AHC Client Version 1.10.00 for Ipnux and Windows! 

O SANC CEientv. 1.10 (1.41 Mb zip-fiiej 
Requirements: 

• 5 MB free disit space on hard drive 

• Tile Java Runtime Environment (JRE) version 5 or iater 

• Linux users require bash 3.0 or iater 

instructions: 

• Ail needed instructions are inciuded in downioad zip-f[ie -- see README.txt 

• For more infbrmation see Users Support. 



Figure 1: SANC Project Download 



2.2 SANC tree for Drell— Yan processes 

The desired stage of development of the SANC tree for any process is the tree with four FORM modules at 
the tip of the processes branches. In Fig. [2] we show the location of the 2/ 2/ CC and NC sub-processes 
at the SANC tree. 

Moving along the menu sequence SANC EW Processes 4 legs ^ 4f ^ Neutral Current 
^fl fl ^ f f (FF,HA,BR,MC) the user can see a hst of the Standard SANC FORM Modules , the 
scalar Form Factors (FF); then at the second module, Helicity Amplitudes (HA); then at the third 
module, the integrated bremsstrahlung (BR); and finally at the fourth module, the fully differential 
bremsstrahlung (MC). 

The same way exists for the branch of Charged Current ^fl fl' ^ f f (FF,HA,BR,MC) 
standard SANC modules, see Fig. [2] 

The FORM modules compute online the FF, HA, BR and MC contributions of the corresponding 
partonic sub-processes. 

Each of these modules in turn can be opened, compiled and run as explained in the User Guide_v 1.10 
(see project homepage). 

^ To install and run SANC client one should have the Java Runtime Environment (JRE) at least version 5.0 Update 5 
installed, see section Minimum System Requirements of the Download page at the SANC project homepage. 
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Figure 2: SANC tree for the Drell-Yan processes: Charged and Neutral Current case 



For more details of running the existing accessible menu, one can refer to the description of the SANC 
system in Ref. [14] . 

Let us consider the types of modules. Each FORM module produces an output, which is in turn an 
input for creation of the corresponding FORTRAN code by the software s2n. 

• FF and HA 

In Ref. [HI we presented the Covariant (CA) and Helicity Amplitudes (HA) for fififf — > NC and 
CC processes, with all 4-momenta being incoming for any of its cross channels s, t or u. The expressions 
for the CA and HA (see Eq. (30) and (33) of the last reference) of these processes are written in terms of 
scalar FF's. 

• BR and MC 

The BR module computes the soft and inclusive hard real photon emission: 

q + q^ £ + £ + J. (3) 

We do not discuss the soft photon contribution here, referring the reader to the system itself. As far 
as hard photons are concerned, we realized two possibilities of the integration over their phase space: 
the semi-analytical one (module BR) and the one by means of a Monte Carlo integrator or generator 
(module MC). 

The MC module provides fully differential hard bremsstrahlung contribution to the partonic cross 
section. The contribution is given in a form suitable for further numerical integration or simulation of 
events in a Monte Carlo generator. 
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2.3 Types of SANC Outputs 

As we mentioned in section [TJ SANC produces two types of outputs: SSFM and MC generators at the 
one-loop level. 

• Standard SANC FORTRAN Modules , i.e. the list of packages at the parton level for download is, 
(see Figini) 

- SANC NC vl.20: SSFMs for the DY NC processes: {uu,dd) ^ i^J-H,ee) [52], SSFMs for the 
processes: ee{uu,dd) — > HZ [15] and NC gluon-induced processes [33] : 



- SANC CC vl.20: SSFM for the DY CC processes: (uu, dd) {fii'^,ei^e) p4] and CC gluon-induced 




One has to emphasize that the CC vl.lO contains first SANC modules for the calculation of NLO QCD 
corrections [3S]. Some results about interplay of QCD/EW correction were reported in talks to ATLAS 
MC group [IH],[1H1,[1S1- 

SANC Modules. 

SANC NC package (Neutral Current processes modules). 

The previous versions of SANC NC pacl;age you can find in SANC Archives. 

O lS/12/2008 SANC NC vl.20 pacl;age (3S4 Kb tgz-file) [last stable version] 
In vl.20 package the modules for new processes are added. Details are in the file CHANGES. 

O 13/02/2008 SANC NC vl.lO package (300 Kb tgz-file) [stable version] 

In vl, 10 package the flflHZ processes are added, For description of the flflHZ processes please refer to the paper D, Bardin 
et a!., Comput. Phys. Comm. 177 (2007) 73S-7S6. Details are in the file CHANGES. 

This package is intended for calculation of the 1-loop radiative correction to Drell-Van Neutral Current processes at partonic 
level, For a Technical Description of this module please refer to the paper Eur, Phys, J. CS4- (200S) 4S1, 

SANC CC package (Charged Current processes modules). 

The previous versions of SANC CC package you can find in SANC Archives, 

O OS/12/2008 SANC CC vl.20 package (160 Kb tgz-file) [last stable version] 
In vl.20 package the modules for new processes are added. Details are in the file CHANGES. 

O 28/10/2008 SANC CC vl.ll package (131 Kb tgz-file) [stable version] 
In vl.ll package some bugs in the QCD soft-virtual part are fiited. Details are in the file CHANGES. 

This package is intended for calculation of the 1-loop radiative correction to Drell-¥an Charged Current processes at partonic 
level, see A, Arbuzov et ak, Eur, Phys, J. C46 (2006) 407, 




processes 




Figure 3: Available versions of the SSFM 
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• Standalone MC generators: 

The sanc_**_foam_vl.** packages are intended for generation of unweighted events of the DY processes 
in NC and CC sector at the hadronic level taking into account the one-loop EW radiative correction 
based on FOAM algorithm 16 , (see Fig. [6]). These generators use the standard SANC Fortran modules 
for calculation of NLO EW corrections as well as LoopTools-2 . 1 [5S] and SancLib-vl-00 libraries for 
evaluation of scalar and tensor one-loop integrals. Also you need ROOT package to be installed at your 
computer to use these generators. 

In the present version of packages we include the possibility to write the output in data files containing 
the event information in the standard Les Houches Accord format [37 in order to organize the transfer 
of information between SANC generators and general purpose programs such as PYTHIA 38J and HERWIG 

m- 

We advice to read INSTALL for installation instructions, UserGuide.txt for a Technical Description 
and readmc-foam for FOAM using. 

Some results obtained with these generators were presented in two talks to ATLAS MC group 17J, 
|18) . Examples of distributions produced with help of MC generator for single-W production are shown 
in Fig. [1 and Fig. El 
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Figure 4: The EW NLO distributions of AfjK and the relative corrections 



Distribution of transverse momentum | 



LO 
NLO 



^0 



45 50 
P^Cii*) [GeV] 




Figure 5: The EW NLO distributions of pj. and the relative corrections 
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SANC Generators. 

Drell-Yan Neutral Current processes generator. 

23/05/2008 SANC NC DY FOAM vLlO package (345 Kb tgz-file) [last stable version] 
In vl.lO package the Les Houches format event files is introduced. Details are in the file CHANGES. 

22/02/2008 SANC NC DV FOAM vl.Ol package {584 Kb tgz-file) [stable version] 
In vl.Ol package the modules are taken from the 5anc_nc_vl.01 package. Details are in the file CHANGES. 

25/12/2007 SANC NC DY FOAM Vi.OO package (583 Kb tgz-flle) 

This SANC generator is intended for generation of unweighted events of the Drell-Yan Neutral Current processes at hadronic level taking into 
account the l-|oop EW radiative correction using the FOAM algorithm. Please read the instructions (in tgz file] for using this package. 

Drell-Yan Charge Current processes generator. 

15/04/2008 SANC CC DY FOAM vl.OO package (235 Kb tgz-fiie) [last stable version] 

This SANC generator is intended for generation of unweighted events of the Drell-Yan Charge Current processes at hadronic level taking into 
account the 1-loop EW radiative correction using the FOAM algorithm. Please read the instructions (in tgz file) for using this package. In i/l.OO 
package the modules are taken from the sanc_cc_vl, 02 package. 

Figure 6: SANC IVIC Generators 

3 Example of the implementation of the SSFM from SANC to 
WINHAC. 

The result of the implementation of SSFM for EW RC at one- loop level is presented in [26] . 

For the description of WINHAC and SANC we refer the reader to the literature: for WINHAC to [40], [41] 
and for SANC to 14J. The goals of the work [IB] were: 

• to check the implementation of SSFM EW NLO modules into the framework of WINHAC Monte 
Carlo event generator; 

• to perform a tuned comparison of two codes: 

1. the standard SANC integrator DY_CC_VEGAS with a modified treatment of ISR QED corrections; 

2. the modified WINHAC, upgraded with the SANC electroweak modules and downgraded to the 0{a) 
QED corrections. 

We reached the agreement between the WINHAC MC event generator and the SANC MC integrator for 
the 0{a) EW corrections to the CC Drell-Yan process at the level of ~ 0.025%. 

As an example we demonstrate the comparison of the transverse mass distributions for the muon 
channel in two popular EW parametrization schemes (a and Gfj,, see e.g. Ref. [6]) with "simplified bare" 
cuts, see Fig. [7] The two upper figures show the RC quantity Sew in %, while the two lower figures show 
absolute deviations W— S between the two calculations, for details see Ref. [26| . 
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Figure 7: The EW NLO distributions of for the muon channel in two schemes and the absolute 
deviations between WINHAC and SANC, A W — S, both quantities in [%]. 

4 User guide 

Here we present the technical description of the SANC NC and SANC CC packages — v. 1.20, intended 
for calculation of the total and differential da /d cos i) cross sections at the partonic level. Conventions, 
enumerations and descriptions of the options of the main flags are also given. 
The packages can be accessed from project homepages [21]. 



4.1 Naming conventions 

In SANC we use naming conventions for fields (or particles) shown in Table [T] where N is the field index 
(here we present only physical particles, omitting ghost fields) and in the columns headed "name" we 
show the names used internally in SANC. All associated parameter symbols are derived from these names. 
Thus the mass, charge and weak isospin of the electron are denoted mel, qel and i3el, respectively, also 
the vector and axial vector coupling constants (vel, ael). 



4.2 Overview of the packages 

As an example, we present an overview of the most recent package sanc_cc_vl.20[l. It contains SSFM of two 
types: the modules cc_*_xy_zv.f(F) for 1-loop level EW sub- processes, which are produced by the s2n pack- 
age of SANC project and are governed by main file main_cc_lloop_vegas.F written "by hand" ; the modules cc- 

^The structure of the package sanc_nc_vl.20 is similar. 
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Table 1: List of fields. 



bosons 


fermions 


QCD 




1st generation 


2nd generation 


3rd generation 




N 


field 


name 


N 


field 


name 


N 


field 


name 


N 


field 


name 


N 


field 


name 


1 


A 


gm 


11 




en 


15 




mn 


19 


Vr 


tn 


23 


g 


gn 


2 


Z 


z 


12 


e 


el 


16 




mo 


20 


T 


ta 








±3 


M/± 


w 


13 


u 


up 


17 


c 


ch 


21 


t 


tp 








4 


H 


h 


14 


d 


dn 


18 


s 


St 


22 


b 


bt 









qcd_*j<y_zv.f for 1-loop level QCD sub-processes and cc*_23_K__y_z_v.f(F) ( cc*_23_x__y_z_v.f(F)) for gluon- 
induced tree level partonic sub-processes together with corresponding main file main_cc_trees_vegas.F. The 
status of QCD SSF(ORM)M is somewhat different from EW ones; the exist, however, at present time 
are not yet put into the system. Nevertheless, SSF(ORTRAN)M were produced from FORM log-files by 
mean of the standard SANC s2n software. Next, the package contains declaration files — *.h and libraries 
— *.a. The folder SancLib_vl_02 contains the source files comprising the library for various functions, 
including Vegas integration and the library itself libSancLib_vl_02.a. Here and below "x,y,z" and "v" 
stand for the standard SANC field indices, i.e. 12 — e, 13 — u, 14 — d, 16 — fi, 23 — g, etc. In Table [2] 
we summarize available in sanc_cc_vl.20 1-loop level and in Table [3] all gluon-induced tree level partonic 
sub-processes, renumbered by the Process IDentifier (PID): 

Table 2: Process IDentifier (PID) and available at 1-loop level CC DY sub-processes. 



PID 


xy_zv 


Sub- 


process 


1 


1314_1112 


u + d - 


De + e~ 


2 


1413.1211 


d + u - 


e+ + 1^^ 


3 


1314_1516 


u + d - 


+ ^1- 


4 


1413.1615 


d + u - 


fi+ + 1^^ 


5 


1314_1920 


u + d - 


Dr + 


6 


1413_2019 


d + u - 





The total set of files inside the package is: 

Instruction files : 
README 

RELEASE-NOTES 

CHANGES 

LICENSE.TXT 

INSTALL 

FILES 

Declaration files : 
s2n_declare . h 
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Table 3: Process IDentifier (PID) and available at tree level gluon-induced sub-processes. 



PID 


23_x__y_z_v 




Sub 


-process 




PID 


(x_23__y_z_v) 




Sub 


-process 




7 


23.13._14_11.12 


9~ 


h-u ^ 


d- 




e+ 


8 


14_23__13_12_11 


9- 


^d^ 


u - 


'r Ve + 


e+ 


9 


23.14..13.12.11 


9- 


^d^ 


u - 


l-e- H 


- i^e 


10 


13.23..14.11.12 


9- 


'rU^ 


d- 


l-e- H 


- Ve 


11 


23.13..14.15.16 


5^ 


- M — > 


d^ 


- -f- 




12 


14.23..13.16.15 


5^ 


-d^ 




- i^A* + 




13 


23.14..13.16.15 


5^ 


-d^ 


u-\ 


-^JL- H 




14 


13.23..14.15.16 


^?^ 


- U —t 


d^ 


-/i- H 


\- Vf_, 


15 


23.13..14.15.16 




- U ^ 


dH 




T+ 


16 


14.23..13.16.15 




-d^ 


u-\ 




T+ 


17 


23.14..13.16.15 


.9H 


-d^ 


u-\ 


-T- H 


- Vr 


18 


13.23..14.15.16 


.9H 


- u ^ 






- Vr 



Initialization and various input files: 

s2n_init . f 
sanc_input . h 
leshw_input . h 
tev41hcw_input . h 
leshw2007_input . h 

SSFM originated from 

cc_f f _xy_zv.F (FF) 
cc_si_xy_zv. f (HA) 
cc_ [qcd_] br_xy_zv . f (BR) 

[this file contains three SSFM (subroutines) for EW case and two for QCD case 

cc_bo_xy_zv( . . . ) , cc_br_xy_zv( . . . ) , cc_ha_xy_zv_lspr ( . . . ) 
cc_qcd_br_xy_zv( . . . ) , cc_qcd_ha_xy_zv_lspr ( . . . )] 
cc_ [qcd_]ha_xy_zv.f (MC) 

Main file: main.cc.lloop.vegas.F 

As a rule of the SANC approach, we subdivide the EW RC into the virtual (loop) ones, the ones due to 
soft photon emission, and the ones due to hard photon emission with the aid of the soft-hard separator 
Lo. For all SANC processes we demonstrate the numerical independence off this auxiliary parameter. The 
adopted form of presentation of the differential cross section at the one-loop level in obvious notation is: 

^^l-loop ^ ^^Born ^ ^^Subt ^ rf^Virt+Soft ( - ^ ^ dG^''"^(G)). (4) 

The second term stands for subtraction of coUinear quark mass singularities. It may be or 
^^MS(Dis)^ correspondingly (see [26] and [34]). At the partonic level only MS option is realized: 

dc^ ^db- Aa^, (5) 

where 

Aa^ - lim [ AaV'^'+s°ft {c) + Aa^^''^ (w)] ^^ (6) 
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It is described in detail in Refs. [31], [32] ■ At the hadronic level, where Eq. ^ is convoluted with 
PDF's in the usual way, the other two "subtraction" options (YFS or DIS) were used. The terms, made 
up of the virtual (loop) ones — da^"^ , the ones due to soft photon emission dcr®°^* , and the ones due 
to hard photon emission da^^'^'^ are subdivided into ISR, IFI and FSR in accordance with VF-splitting 
techniques [B]. An auxiliary parameter uj separates the soft and hard photonic contributions. 

The steps of calculation in the main_cc_lloop_vegas.F file are in accordance with Eqs. 

• step of declaration and initialization , followed by a call Processlnit (pid), see Section l473l for descrip- 
tion of various options of flags; 

• step born is realized by flag iborn=l, 

^g.Born jg computed by integration over c, see Ref. [32] , of the function_ew_lc , 
via call ProcessBorn (...,born,...) 

— and SSFM call cc_br_xy_zv(...,born,...) 

• step virt+soft and all subsequent ones are realized by flag iborn^O, 
^^Virt+Soft jg computed by integration over c of the function_(ew,qcd)_lc, 
via call ProcessVirt_(EW,QCD) (...,virt,...) 

— and Virt(Weak) by SSFM call cc_si_xy_zv(...,sigma) , 

(inside this module there exists a call to SSFM cc_ff_xy_zv(-s,-t,-u)) 

— and Virt+Soft(qED,qCD) 

* for QED processes by SSFM call cc_br_xy_zv(...,soft,...) 

* for QCD processes by SSFM call cc_qcd_br_xy_zv(...,soft,...) 

• step brdq, i.e. (Act)^^, 

via call ProcessBrdq_(EW,QCD) (...,brdq,...) 

a) first part, i.e. (A^virt+Soft)MS ^ 

— is computed by SSFM call cc_bo_xy_zv(...,born), multiphed by ((JSV^ms — (jq analog of 
Eq.(19) of [32] 

b) second part, i.e. {da^'^'^'^)^^ is computed by integration over ^ of the function_(ew,qcd)_ls 

— and by SSFM call cc_bojcy_zv(...,bornk), using Eq.(12) of [34] 

• step hard , i.e. da^'^'^'^ / ds' is computed by integration over s' of the function_(ew,qcd)_lspr 
via call ProcessHard_(EW,QCD)_lspr (....hard,...) 

— and SSFM call cc_[qcd_]ha_xy_zv_lspr (...,hard) 
or alternatively 

• step hard , i.e. da^^^'^ / d^^"^^ is computed by integration over 4d-phase space of Eq. (6)-(7) of [32 
of the function_(ew,qcd)_4d 

via call ProcessHard_(EW,QCD) (....hard,...) 

— and SSFM call cc_[qcd_]ha_xy_zv(...,hard) 
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Main file: main_cc_trees_vegas.F 

The presentation of the differential cross section for this case is trivial: 

Here the Born-term describes the tree level cross section of one of processes of Table [3] The subtraction 
of collinear final quark mass singularities for the time being is applied in integrators and generators and 
is not implemented at the parton level. 

4.3 Flag options 

pid(I) — choice of the 1-loop level sub-process, see Table El 

1=1,.. .,6 
or of the tree level ones of Table [31 
pid(I) 1=7,. ..18 

iqed(I) — choice of calculations for QED correction: 

1=0 without QED correction 

1=1 with all QED correction 

1=2 with ISR QED correction 

1=3 with IFI QED correction 

1=4 with FSR QED correction 

1=5 with IFI and FSR QED correction 
iew(I) — choice of calculations for EW correction: 

1=0 without EW correction 

1=1 with EW correction 
iqcd(I) — choice of calculations for EW correction: 

1=0 without QCD correction 

1=1 with QCD correction 
iborn(I) — choice of scheme of calculations of the partonic cross section: 

1=1 only Born level 

1=0 Born + 1-loop virtual corrections 
gfscheme(I) — choice of the EW scheme: 

1=0 Q;(0)scheme 

1=1 Gf scheme 

1=2 Gp scheme 
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ilin(I) — choice of the Unearization at the calculation of the partonic cross section: 
1=0 without linearization 

1=1 with linearization, i.e. neglecting spurious terms 0{a'^) 
ifgg(I) — choice of calculations of photonic vacuum polarization J^gg-. 
I=-l 
I== 1 
1= 1—1 + kJCgg 
1= 2 - 1/(1 - k^gg) 

2 

with k = 7^— T . 

lOTT 

ihard(I) — types of the hard bremsstrahlung phase-space integrations: 

1=1 integration over s' 

1=4 4d integration 
isetup(I) — choice of the setup: 

1=0 Standard SANC 

1=1 Les Houches Workshop, 2005 

1=2 TcV4LHC Workshop, 2006 

1=3 Les Houches Workshop, 2007 

5 Conclusions 

In this paper we have described the Standard SANC FORTRAN Modules and presented examples of their 
application: 1) they were used in the packages at the parton level for quick studies of different features of 
the given sub-proccsscs: estimates of effects due to variations of input parameters, elcctrowcak schemes, 
interplay of different RC contributions (EW-QCD, initial and final state radiation) etc.; 2) they were 
implemented in the MC generator WINHAC, increasing thereby its potential possibilities; 3) finally, they 
were used in the SANC Monte Carlo integrators and event generators that provide distributions of the 
final state particles with full kinematics, which was interfaced with parton showering codes (PYTHIA and 
HERWIG) and the events can be further processed through the whole experiment simulation environment. 
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